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Abstract 

We study scalar (qq) and (qq) 2 correlation functions in the instanton liquid model. We show that 
the instanton liquid supports a light scalar-isoscalar (sigma) meson, and that this state is strongly 
coupled to both (qq) and (qq) 2 - The scalar-isovector ao meson, on the other hand, is heavy. We 
also show that these properties are specific to QCD with three colors. In the large N c limit the 
scalar-isoscalar meson is not light, and it is mainly coupled to (qq). 
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I. INTRODUCTION 



There is a longstanding controversy concerning the structure of scalar-isoscalar mesons 
in QCD . One view is that scalar mesons are p-wave qq bound states 

with masses in the (1.0-1.5) GeV range. According to this view the large effects seen in the 
I = tttt and 1 = 1/2 ttK channel are due to t-channel meson exchanges and other coupled 
channel effects. The competing view is that there are genuine scalar resonances below 1 
GeV, and that the structure of these states is quite different from ordinary quark model qq 
states. 

Jaffe suggested that the unusual properties of the light scalar mesons could be explained 
by assuming a large (qq) (qq) admixture . He observed that the spectrum of the flavor 
nonet obtained by coupling two anti-triplet scalar diquarks is inverted as compared to a 
standard qq nonet, and contains a light isospin singlet, a strange doublet, and a heavy 
triplet plus singlet with hidden strangeness. This compares very favorably to the observed 
light sigma, the strange kappa, and the heavier a (980) and /o(980). It also explains why 
the ao and fo are strongly coupled to KK and ttt]. Strong correlations in the anti-triplet 

in QCD and have also been invoked in order to 



scalar diquark channel are quite natura 

understaud the structure of baryous QQ aud deuse uaryou.e matter QQ 

In this work we study the nature of scalar-isoscalar mesons in the instanton model |14ill£ 
There are several reasons why instantons effects are relevant. Scalar-isoscalar mesons have 
vacuum quantum numbers and are directly related to fluctuations in the chiral condensate. 
There is an impressive amount of evidence from lattice calculations that instantons are 
responsible for chiral symmetry breaking in QCD 12, 3- These calculations also suggest 
that the approximate zero modes responsible for chiral symmetry breaking are small in size. 
This implies that the chiral condensate is very inhomogeneous and that fluctuations in the 
scalar channel are large. In addition to that, if there is a large (qq)(qq) admixture in the 
wave function of scalar mesons, both flavor mixing and off-diagonal (qq) — (qq) 2 mixing have 
to be large. The instanton induced (qq) Nf interaction provides a natural mechanism for 
these effects. 

In order to study the nature of scalar-isoscalar mesons we compute both (qq) and (qq) (qq) 
correlation functions. We show that there is a large scalar resonance signal in both the 
diagonal (qq) and (qq)(qq) correlators as well as the off-diagonal (qq)-(qq)(qq) correlator. 
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This implies that there is indeed a large four quark admixture in the scalar-isoscalar meson. 
We have to note, however, that in a relativistic quantum field theory this statement is 
somewhat hard to quantify. We can compute the coupling of the (qq) and (qq) (qq) operators 
to the scalar resonance, but these coupling constant have different dimensions and it is not 
clear how to compare them. 

We address this problem by studying scalar correlators for different numbers of colors. 
We expect that in the large N c limit the scalar meson becomes a "normal" quark model state 
which is dominated by the (qq) component. This implies that there is no scalar resonance 
signal in the four quark channel, and that the scalar-isoscalar and scalar-isovector correlators 
become indistinguishable. We show that this tendency is clearly seen in ratios of correlation 
functions computed in QCD with different numbers of colors. 

This paper is organized as follows. In Sect. [H] we define (qq) and (qq)(qq) scalar cor- 
relation functions. In Sect. II I II we discuss the short distance behavior of these correlators 
and investigate the importance of unphysical contributions in the quenched approximation. 
In Sect. II VI we present numerical results from quenched and unquenched simulations of the 
instanton liquid model. 

II. SCALAR CORRELATORS 

The lowest dimension operator with the quantum numbers of the sigma meson is the 
scalar (qq) operator 

3(qq) = -^(uu + dd). (1) 
The corresponding correlation function is (see Fig. [IJ 

U {qq) (x) = 0(x) - 2T(x) = (Tr [S(0, x)S(x, 0)]> - 2(Tr [S(0, 0)] Tr [S(x, x)}), (2) 

where (.) is an average over all gluonic field configurations, S(0, x) is the full quark propaga- 
tor in a given background gluonic field, and the trace is taken over Dirac and color indices. 
0(x) and T(x) denote the one-loop (connected) and two-loop (disconnected) terms in the 
correlation function. Note that the quark line connected term is equal to the correlator of 
the scalar-isovector current, U a Jx) = O(x). Lattice calculations of the correlation function 
equ. (0) have been reported in 
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There are several ways to construct scalar-isoscalar four quark operators. One possibility 
is to couple two quark-anti-quark operators. We expect the most attractive channel to arise 
from the coupling of two color-singlet, flavor-triplet pseudoscalar operators, 



3{qqy = {qi^ a q){qi^ a q). 



(3) 



The corresponding correlation function is equal to the 1 = two-pion correlator 21 

n /0 (x) 

where 



(4) 



D(x) 
C{x) 
A(x) 
G{x) 



((Tr[5(0,x) 75 5(x,0)75]) 2 ) 
(Tr([5(0,a;)755(x,0)75] 2 )) 
(Tr [S(0, x) l5 S(x, x) l5 S(x, 0) 75 5(0, 0) 75 ]) 



(Tr (5(0,0)75) 



Tr 



(S(x,x)^ 5 ) 2 ) 



(5) 
(6) 
(7) 
(8) 



are the four quark correlators shown in Fig. El The 1 = 2 two-pion correlator is given by 



U I2 (x) = D(x) - C(x) 



(9) 



and does not involve any disconnected diagrams. The mixed {qq)-{qq) 2 correlation function 
is given by (see Fig. |HJ) 



where 



U^gy I0 (x) = 7y (j( qq -)(0)j(qq)s(x)) = (A mix (x) - G mix {x)) 



A mix {x) = (Tr [S(0,x)-f 5 S(x,x)'y 5 S(x, 0)]) 



G mix {x) = (Tr[S'(0,0)]Tr (S(x,x) l5 ) 2 > 



(10) 



(11) 
(12) 



Instead of coupling two quark-anti-quark operators we can construct scalar-isoscalar four 
quark operators by coupling a diquark and an anti-diquark operator. The most attractive 
channel is expected to arise from two color anti-triplet and isospin singlet scalar diquarks 



JW = e abc e ade (q b C 75 r 2 q c )(q d C 75 r 2 q e ) 



(13) 
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This operator can be Fierz-rearranged into a (qq) 2 operator. The result is 

j(dq)' = \{{qi^ a q)(q^f a q) + {qf a q){qf a q) + (q-f^f a q)(q-f^f a q) 

+ (qi5l„f a q) (q-f 5 l„f a q) - - {qa^f a q) {qa^f a q) } ( 14) 

where f a = (r,i). Note that the flavor structure of equ. (fT4*|) is that of a determinant, 
{qf a q){qf a q) = (uu)(dd) — (ud)(du). Equation (JTljl contains both attractive channels, such 
as 7T7T and /o/o, as well as repulsive channels, such as rj'r]' and a^aQ. In the following we will 
concentrate on the tttt channel equ. (J3J). 



III. SCALAR CORRELATORS IN THE INSTANTON MODEL 

In this section we wish to study instanton contributions to the scalar correlation functions 
equ. (|2|4|9jl and equ. pOjl . At short distances only the effect of the closest instanton has to 
be taken into account. The main instanton contribution is related to the zero mode term in 
the fermion propagator 

S(x,y)^ Mx - z)4{y - Z \ (15) 
m* 

where i[)o(x) is the zero mode wave function and z is the location of the instanton. For an 
isolated instanton m* = m q where m q is the current quark mass. In an ensemble of instantons 
m* is an effective mass m* ~ TTp(2/N c ) 1 ^ 2 (N/V) 1 ^' 2 . Here, p is the average instanton size 
and N/V is the average instanton density. Instanton zero modes make equal contributions 
to the one and two-loop terms in the scalar correlation function equ. (J2J), Oi nst = T inst < 0. 
As a consequence the one-instanton contribution is attractive in the a channel and repulsive 
in the ao channel. Averaging over the position of the instanton we get 

n»^/*»W^| F g( I $ + ^ 8 ^)j, (16) 

where n(p) is the instanton size distribution and £ 2 = x 2 /(x 2 + 4p 2 ). The one-instanton 
contribution can be resummed using the random phase approximation (RPA) 24 , 2^| |^ 



The loop and vertex functions C s and T s are given by 

CM) - 4iV c [-) J (27r)4 { M 2 +P 2 )(M 2 +p 2 ) ' (18) 
rn f d 4 p (M.M^jM^-p.-p,) 

s{q> ' J (2vr) 4 (M 2 + p 2 ){M 2 + p 2 ) ' 1 ' 



where pi = p + q/2, pi = p — q/2 and Mi j2 = M(pi i2 ) are momentum dependent effective 
quark masses. In the instanton liquid model the momentum dependence of the effective 
quark mass is governed by the Fourier transform of the instanton zero mode profile. The 
effective quark mass at zero virtuality, M(0), is determined by a Dyson-Schwinger equation. 
For typical values of the parameters p ~ 0.3 fm and (N/V) = 1 fm -4 we find M(0) ~ 350 
MeV. 

Scalar correlation functions in the RPA approximation are shown in Fig. (j3J). We observe 
that the isoscalar and isovector correlation functions are indeed very different. We also note 
that the isovector correlator is unphysical for x > 0.7 fm. This is an artifact of the quenched 
approximation and can be understood in terms of quenched irr]' intermediate states. The 
quenched a correlation function receives unphysical tj't]' contributions, but these terms are 
not included in the random phase approximation. The scalar-isoscalar (cr) correlator shows 
a clear resonance signal. A simple pole fit gives m a = 540 MeV. We observe that the a and 
do correlation functions become more similar as the number of colors increases. However, 
in the RPA approximation the difference remains finite and large in the limit N c — > oo. As 
explained in |22| this is related to a failure of the RPA approximation when applied to the 
large N c instanton liquid. 

Let us now study instanton contributions to the two-pion correlation functions equ. (J3J 
and (JOJ). The first possibility is that all four quarks propagate in zero mode states. In QCD 
with Nf = 2,3 flavors this contribution violates the Pauli principle and has to vanish. This 
is reflected by the fact that in both the 1 = and 1 = 2 correlation functions the sum of 
the coefficients of the four different contractions, D,C,A and G, vanishes. We note that 
for light quark masses it is dangerous to drop the disconnected contributions to the 1 = 
correlator, an approximation sometimes referred to as the quenched valence approximation. 
In particular, there is a large positive zero mode contribution to the connected 1 = 
correlation function which is not present in the full 1 = correlator. 

Contributions with an odd number of fermion zero modes violate chirality. This leaves 
terms with two zero mode propagators, see Fig. |H1 These terms only appear in the direct 
and crossed diagrams. The one-instanton contribution to the direct diagram is D SIA (x) = 
211° (x)U^ IA (x), where 11° (x) = 3/(7r 4 x 6 ) is the free quark contribution and U^ IA (x) = 
(x) is the single instanton contribution to the pion correlation function. It is clear 
that D SIA (x) simply corresponds to two non-interacting pions. Interactions arise from the 
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crossed term. We find 

IO*) = 2IT°(*)ie 4 (x) + in°(x)nf^(x) (20) 
US A (x) = 2Ul(x)U s n IA (x) - ll^{x)I^ IA {x). (21) 

This implies that the one-instanton term is attractive in the 1 = channel and repulsive in 
the 1 = 2 channel. One can try to resum this interaction in order to study the irir interaction 
at intermediate and lon g di stances. For the one-gluon exchange interaction, this problem 
was recently studied in 28[. At very low momenta, the nir interaction is constrained by 
chiral symmetry. Weinberg predicted the s-wave scattering lengths in the 1 = and 1 = 2 
channels j^J 

_ 7m w 2 

°"32vr/r °" 167T/2- {U) 

We emphasize, however, that our objective in this work is not a determination of the scat- 
tering lengths, but a study of scalar resonances. In Sect. II VI we will describe how to sum 
the instanton induced interaction to all orders by performing numerical simulations of the 
I = 0, 2 7T7T correlation functions in the instanton model. 

At the one-instanton level there is no contribution to the the mixed (qq)-{qq) 2 correlation 
function. This is a consequence of the fact that in Nf = 2 QCD not all quarks can propa- 
gate in zero mode states. In QCD with three flavors there is a flavor mixing (uu)(dd)(ss) 
contribution. 

Finally, we have to understand the effect of the quenched approximation on the isoscalar 
two-pion correlation function 30, see Fig.[7| The main artifact in the quenched approx- 
imation arises from unphysical rj'rf intermediate states. The long distance, low momentum 
part of the quenched rf propagator is of the form 

Mff) = tx-t ~ ^r^ m l^r^> ( 23 ) 

where the second term is the so-called double-hairpin contribution, Wq ~ mi, is related to 
the rj' mass and A,r is the coupling of the pion to the pseudoscalar current. Note that the 
double hairpin term is negative in euclidean space. Using equ. (}2"3*j) we see that the rj'r)' 
contribution to the I = tttt correlation function contains a one double-hairpin term which 
is negative and a two double-hairpin term which is positive in euclidean space. Numerical 
results for the coordinate space correlation function are shown in Fig. |SJ The correlation 
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function is normalized to the square of the pion correlator. For comparison we also show the 
contribution of a scalar o resonance with mass m a = 500 MeV and coupling constant X a = 
(350 MeV) 5 . We observe that the quenched approximation leads to unphysical contributions 
that are repulsive at intermediate distances x ~ 0.5 fm and large and attractive at large 
distance. We also observe that these effects will mask the presence of a sigma resonance, 
unless the sigma coupling is very large. 



IV. NUMERICAL RESULTS 



In this section we present numerical results obtained in quenched and unquenched simu- 
lations of the instanton liquid model The main assumption of the instanton liquid 
model is that the QCD partition function 

N f 

Z = I DA^ exp(-S) JJdet(i^ + im f ) (24) 

/ 

is dominated by classical gauge configurations called instantons. Instantons in QCD with 

N c colors are characterized by 4N C collective coordinates, position (4), size (1), and color 

orientation (4iV c — 5). An ensemble of instantons is governed by the partition function 

i , n t +n a N f 
Z =H WW II [dM(pt)]exp{-S)lldeb{iI} + im f ), (25) 

where Q = (zi,pi,Ui) labels the collective coordinates and d(p) is the instanton measure. 
Correlation functions are computed from the quark propagator in a given instanton config- 
uration. The propagator is determined by numerically inverting the Dirac operator in the 
space of approximate zero modes. The short distance part is added perturbatively. 

Most of the results presented below were obtained from simulations at a quark mass m q = 
20 MeV in a euclidean box with dimensions (2.83 fm) 3 x 5.66 fm. We have not systematically 
studied finite volume corrections. Finite volume effects on the two-pion correlator can be 
used in order to determine pion scattering lengths 34]. Our aim in this work is not a 
determination of the tttt scattering length, but a study of tttt resonances. 

Scalar-isoscalar and scalar-isovector correlation functions are shown in Fig. We note 
that the scalar-isoscalar a and scalar-isovector correlators behave very differently. The 
ao correlator is very repulsive. In the quenched approximation the correlation function even 
becomes negative. This behavior can be understood in terms of quenched r/V intermediate 
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states. In full QCD the correlator is still quite repulsive. We extract a resonance mass 
m ao ~ 1 GeV. The sigma correlator, on the other hand, is very attractive, both in quenched 
and full QCD. Because of the large vacuum contribution ((qq) 2 ) it is hard to extract the 
mass. Our results are consistent with a light scalar m a = (0.55 — 0.70) GeV. 

In Fig. El w e show / = 0, 2 tttt correlation functions in the quenched approximation. 
The results are normalized to the square of the pion correlation function. We observe 
that the I = 2 tttt correlator is indeed repulsive, the correlation function is suppressed as 
compared to the square of the pion correlator. The connected part of the 1 = correlator 
is very strongly enhanced and has an unusual shape. This is due to the effect discussed 
in Sect. 11111 The connected 1 = correlator receives a large Pauli principle violating zero 
mode contribution. The full 1 = correlator has a vacuum contribution ((qr a q) 2 ) 2 that 
needs to be subtracted. The subtraction involves large cancellations and the result has large 
statistical uncertainties for x > 0.5 fm. We note that the subtracted correlation function is 
repulsive at intermediate distances x ~ 0.5 fm. As explained in Sect. Mil this is an artifact 
of the quenched approximation. 

Unquenched tttt correlation functions are shown in Fig. ^2 We observe that the 1 = 2 
correlation function is not strongly modified as compared to the unquenched calculation. 
The same is true for the connected 1 = correlator. However, the full (subtracted) 1 = 
correlator is strongly modified. The correlation function is now attractive at intermediate 
distances x ~ 0.5 fm. The data are not sufficiently accurate to determine the mass of a 
resonance, but the resonance parameters can be estimated if the data are combined with 
the diagonal (qq) and off-diagonal (qq) — (qq) 2 correlators. We find m a = 600 MeV and 
K = (330 MeV) 5 . 

In Figs. n"2"imi we show the behavior of the correlation functions for different numbers of 
colors N c = 3, . . . , 6 We find that the difference between the scalar-isoscalar and scalar- 
isovector correlation functions disappears as the number of colors increases. For N c = 6 the 
mass of the sigma is m CT > 1 GeV. The physical reason for the change in the sigma meson 
correlation function is the fact that the chiral condensate becomes more homogeneous as 
the number of colors increases. The main change in the scalar-isovector correlation function 
is the disappearance of quenched t]'tt intermediate states. We expect that the unquenched 
ao mass is fairly independent of the number of colors. 

Fig-EHshows that the off-diagonal {qq) — {qq) 2 correlation function decreases dramatically 
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as the number of colors increases. The results are consistent with the idea that the off- 
diagonal correlation function vanishes in the large N c limit. The I = 2 tttt correlator is only 
weakly affected by the number of colors. We observe, however, that the 1 = 2 correlation 
function tends towards the non-interacting tttt correlator. The connected 1 = correlator, 
on the other hand, shows a very dramatic decrease towards the non-interaction two pion 
correlation function as the number of colors increases. 



V. SUMMARY 

In this paper we studied scalar (qq) and (qq) 2 correlation functions in the instanton model. 
We showed that the structure of the scalar-isoscalar a and scalar-isovector do mesons is very 
different. The correlation function in the ao channel is very repulsive while the a channel 
is strongly attractive. The scalar (qq) correlator is consistent with a light sigma resonance 
with a mass m CT = (550 - 700]MeV fl, Q, Q, Q- 

A light sigma resonance also appears in 
Nambu-Jona-Lasinio models |3J|. However, the large OZI violation in the scalar channel is 
characteristic of the instanton induced interaction. 

We also showed that the sigma couples strongly to the 1 = (qq) 2 current. Observing this 
effect not only requires a calculation of both the connected and disconnected diagrams, but 
also the inclusion of the fermion determinant. If the calculation is restricted to connected 
diagrams, or performed in the quenched approximation, the scalar resonance is masked by 
large unphysical contributions. 

We showed that the special properties of the sigma meson, its small mass and large 
coupling to (qq) 2 states, are specific to QCD with three colors. If the number of colors is 
increased the sigma becomes an "ordinary" meson with a mass in the 1 GeV range and 
small coupling to (qq) 2 states. The physical reason is that the chiral condensate is very 
inhomogeneous in N c = 3 QCD, but becomes more homogeneous as the number of colors 
increases. The large inhomogeneity of the condensate in N c = 3 QCD reflects the small size 
of the instanton. The chiral condensate becomes more homogeneous as N c grows because 
the instanton liquid is more dense. As emphasized in (^tJ this does not necessarily imply 
that instantons overlap strongly, since the instanton density grows like the volume of the 
color group. 

There are many interesting problems that remain to be studied. In this paper we concen- 
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trated on scalar mesons in Nf = 2 flavor QCD. One of the remarkable aspects of {qq){qq) 
states in Nf = 3 QCD is their flavor structure. In particular, one would like to verify 
that there is hidden strangeness in the heavy scalars. In addition to the meson-meson or 
diquark-anti-diquark interaction we would also like to study diquark-diquark correlations. 



t 



re spectrum of 



This problem is relevant to the structure of dense baryonic matter 
{qq) 2 q pentaquark states (37 1, and the possible existence of the H-dibaryon 
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FIG. 1: Quark line diagrams contributing to the scalar (qq) correlation function, (O) "one-loop", 
and (T) "two-loop" . The lines are quark propagators in a gluonic background field. 




FIG. 2: Quark line diagrams contributing to the two-pion correlation function, (D) "direct", (C) 
"crossed", (A) "single annihilation", and (G) "double annihilation" or "glue". 
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FIG. 4: Single instanton contributions to the scalar correlation function. The lines crossing the 
shaded area correspond to zero mode propagators, and L, R label the chirality of the quark. The 
anti-instanton contribution corresponds to L *-+ R. 
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FIG. 5: Scalar-isoscalar (a) and scalar-isovector (ao) correlation functions for different numbers 
of colors N c calculated in the RPA approximation to the instanton liquid model. The correlation 
functions are normalized to free field behavior T1q(x) ~ 1/x 6 . The upper/lower curves correspond 
to the cf/clq, respectively. 
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FIG. 6: Instantons contributions to the nir correlation functions. At the one-instanton level, there 
are no contributions to the single and double annihilation diagrams. 




FIG. 7: Diagrams that lead to quenching artifacts in the two-pion correlation function. The 
diagram on the left and right show the two and one double-hairpin terms. 
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FIG. 8: Quenched r/'rj' contribution to the 1 = tttt correlation function normalized to the square 
of the pion correlation function. The curves labeled (lh) and (2h) show the one and two double- 
hairpin diagram, see Fig. Q For comparison, we also show the contribution of a scalar resonance. 
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FIG. 9: Scalar-soacalar a and scalar-isovector ao correlation functions in the instanton model. 
The correlators are normalized to free field behavior ~R§{x) ~ 1/x 6 . The labels (qu) and (unq) 
refer to the quenched and unquenched correlation functions. 
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FIG. 10: 1 = and 1 = 2 two-pion correlation functions measured in quenched simulations of 
the instanton liquid model. The correlators are normalized to the square of the pion correlation 
function. The curve labeled (conn) shows the connected part of the 1 = correlator, (full) is the 
complete 1 = correlator, and (subtr) is the complete correlator after the vacuum contribution is 
subtracted. 
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FIG. 11: 1 = and 1 = 2 two-pion correlation functions measured in unquenched simulations of 
the instanton liquid model. Curves labeled as in Fig. 
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FIG. 12: Scalar-isoscalar (a) and scalar-isovector (ao) correlation functions for different numbers of 
colors N c . The correlators are normalized to free field behavior, n°(x) ~ N c /x e . The upper/lower 
curves correspond to the cx/ao, respectively. The data shown in this figure were obtained from 
quenched simulations of the instanton liquid model. 
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FIG. 13: Off-diagonal (qq)-(qq) 2 correlation function for different values of the numbers of colors 
N c . The correlation functions are normalized to the free (qq) correlator. The data shown in this 
figure were obtained from quenched simulations of the instanton liquid model. 



23 




FIG. 14: Connected 1 = and 1 = 2 two-pion correlation functions for different values of the 
numbers of colors N c . 
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